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ABSTRACT 

We present here studies of collapse of purely baryonic Population III objects with 
masses ranging from IOM 0 to 10 ®Mq. A spherical Lagrangian hydrodynamic code 
has been written to study the formation and evolution of the primordial clouds, from 
the beginning of the recombination era (zr-ec ~ 1500) until the redshift when the 
collapse occurs. All the relevant processes are included in the calculations, as well 
as, the expansion of the Universe. As initial condition we take different values for 
the Hubble constant and for the baryonic density parameter (considering however a 
purely baryonic Universe), as well as different density perturbation spectra, in order 
to see their influence on the behavior of the Population III objects evolution. We find, 
for example, that the first mass that collapses is 8.5 x IO^Mq for h = 1, U = 0.1 
and 6i = 5p/p = (M/Mo)“^/®(l + Zrec)~^ with the mass scale Mg = 10^®Mq. For 
Mo = 4 X IO^^Mq we obtain 4.4 x 1O'^M0 for the first mass that collapses. The 
cooling-heating and photon drag processes have a key role in the collapse of the clouds 
and in their thermal history. Our results show, for example, that when we disregard 
the Compton cooling-heating, the collapse of the objects with masses > 8.5 x 1O^M0 
occurs earlier. On the other hand, disregarding the photon drag process, the collapse 
occurs at a higher redshift. 
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1 INTRODUCTION 

The formation of galaxies, or pre-galactic objects, in the 
early Universe is one of the most important fields of study 
in cosmology. It is not completely known, however, how the 
various structures observed have been formed. It is not com¬ 
pletely clear how the density perturbations, that are the 
seeds for the formation of the structures of the Universe, 
were produced and how they evolved to produce the struc¬ 
tures that we observe. 

Some authors argued that for scales below O.lMpc the 
purely gravitational simulations are not completely appro¬ 
priate, due to the fact that non-gravitational physical pro¬ 
cesses are preponderant (Cen 1992). Thus, the more con¬ 
venient manner of studying these scales is to use a hydro- 
dynamic code and to include all the physical processes that 
are important during and after the recombination era, start¬ 
ing the calculations from the beginning of the recombination 
era. As it will be shown in the present paper, we make such 
a study. 

Studies including the cooling-heating processes of the 
gas, in particular of pre-galactic composition, have been con¬ 
sidered since the 1960’s (e.g., Matsuda et al 1969, Hutchins 
et al 1976, Calberg 1981, Palla et al 1983 among others). 
These articles consider pressure free gas clouds collapsing 
and analyzed the effects of the physical processes. Recently, 


some authors considered the pressure of the gas among other 
physical processes in the gas dynamics (de Araujo 1990, de 
Araujo & Opher 1988, 1989, 1994, Haiman et al 1995, Thoul 
& Weinberg 1995 among others). 

For example, de Araujo & Opher (1988) studied the 
collapse of primordial clouds using an isothermal density 
perturbation spectrum. They included the internal pressure, 
the influence of the cosmic background radiation and several 
cooling-heating processes; among them are: photon cooling, 
photon drag, photoionization due to the CBR, collisional 
ionization and a set of equations for the creation and de¬ 
struction of H 2 molecules, as well as the expansion of the 
Universe. All these processes are present during and after 
the recombination era, therefore, they are essential to be 
taken into account. 

In particular, de Araujo & Opher (1989) used a com¬ 
plete set of hydrodynamic equations and assumed that the 
density contrast profile is uniform (“top hat” like) through¬ 
out the calculation. Thus, the resulting system of equations 
is dependent only on time. As a consequence, the velocity 
profile is linear in the spherical radial coordinate. This ap¬ 
proximation could be good but a detailed analysis shows 
that the physical parameters can change their values for 
different points inside the cloud, as we are showing in the 
present work. Also, it is not realistic to consider that the 
cloud, as a whole, will collapse; a possible partial collapse is 
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more realistic. 

In the present article we used the same set of hydrody¬ 
namic equations of de Araujo & Opher (1988). We rewrote it 
to use a spherical Lagrangian hydrodynamic code to follow 
the variables of the fluid. We analyzed the internal structure 
of the cloud during the collapse, using the same spectrum 
of density perturbations used by de Araujo & Opher (1988), 
to verify the differences in the behavior of the physical pa¬ 
rameters when we include the internal morphology. 

The spherical Lagrangian hydrodynamic code (written 
by one of us, ODM) divides the cloud in several concentric 
shells. This number of shells is the optimal combination for 
the variable time step and the number of shells that gives the 
shortest processing time. Our results show that 450 shells are 
enough to model satisfactorily all the clouds here studied. 
To verify if the number of shells is an optimal choice we ran 
a couple of models with 4500 shells, and we see that the re¬ 
sults are the same as for the models with 450 shells. We also 
reproduced the results of de Araujo & Opher (1994) with 
good agreement, for the models that have small pressure 
gradients. The shocks that appear are treated with the in¬ 
clusion of the artificial viscosity of Richmeyer and Newmann 
(Richmeyer & Morton 1967). 

In another article, to appear elsewhere (Miranda 1998), 
all the tests and details of the spherical Lagrangian hydro- 
dynamic code used in the present article are given. 

We begin the calculations at Zrec = 1500 (Tree = 
4000A), the beginning of the recombination era, and we fol¬ 
low the evolution of the cloud until its collapse. The initial 
density perturbation can be given in the form of a power 
law spectrum as considered, for example, by Gott & Rees 
(1975) 



with M the mass of the cloud. Mo is the reference mass and 
Zrec is the redshift of the recombination. We use flo = 0.1 
(baryonic density parameter) and H — 100km s~^Mpc~^ 
(the Hubble constant). 

It is worth noting that the above equation gives the 
rms density fluctuations as a function of the mass scale. 
This equation sets the initial density contrast in our top hat 
profile model. 

In many studies, dealing with the evolution and collapse 
of primordial structures (some of them referred to above), it 
is not specified when (at what redshift) the calculations are 
started. It is also considered that the proto-structure, under 
study, is already detached from the rest of the Universe and, 
as a result, it is not taken into account the expansion of the 
Universe. The merit of de Araujo and Opher’s (1988, 1989, 
1994) works was to take into account all the necessary ingre¬ 
dients to study the evolution of density perturbations but, 
by other hand, they unrealisticaly fixed the density profile of 
the perturbations during all the time of the cloud evolution. 

In order to follow the evolution of density perturbations, 
it is necessary to start the calculations at the beginning of 
the recombination era, as mentioned above, including the ex¬ 
pansion of the Universe and all the relevant processes that 
take place during and after the recombination era. Obvi¬ 
ously, the results of the study will depend also on the density 
parameter, on the Hubble constant and on the density per¬ 


turbation spectrum chosen as initial conditions. However, 
the three above ingredients are not completely known and 
thus, in the present article, we are considering different com¬ 
binations of these three ingredients to follow the influence 
caused in the evolution of the clouds under study. 

In §2 we describe the basic equations, in §3 we discuss 
our results and finally in §4 we present our conclusions. 


2 BASIC EQUATIONS 

The hydrodynamic equations that describe the dynamics of 
the primordial clouds are: 


dp 1 d . 2 s o 


( 2 ) 


the mass conservation equation written in spherical coordi¬ 
nates; where p is the density of the cloud, r the radial coordi¬ 
nate and V the velocity. In our code the continuity equation 
is calculated directly from the grid. The specific volumes of 
the Lagrangian shells give the density of the cloud. 

The equation of motion is given by: 


+ i VP-h V<(.-h ^21^^ 

Dt p rripC 


\ m . 

R{t) 


= 0 


(3) 


where the velocity is u = with the peculiar veloc¬ 

ity of the cloud and — Hr the Hubble flow; the pressure 
of the gas is P = kN pTm.{l + x), with k the Boltzmann con¬ 
stant, N is the Avogadro’s number, Tm. is the temperature 
of the matter inside the cloud, and x the degree of ion¬ 
ization. The gravitational potential is cp, the Thomson 
cross section, b = do-g^/c (with the Stefan-Boltzmann 
constant), Tr is the temperature of the radiation, is 
the proton mass, c the velocity of light, H{t) = R{t)/R{t) 
with R{i) the scale factor and R{t) its time derivative, and 
D/Dt = dt + vdr is the total derivative. 

The energy equation is written as: 


DE _ P Dp 


(4) 


where E is the thermal energy per gram and L the cooling 
function. 

The cooling-heating processes are included in the cool¬ 
ing function L which is given by the summation of four mech¬ 
anisms: 


L — L^-\- -h + La (5) 

The heat conduction is not taken into account due 
to the fact that the clouds, that we are studying, never 
reach temperatures in which heat conduction is important. 
Bremsstrahlung was also not included because the temper¬ 
atures of the matter inside the clouds never reach values 
above of 10"^K. 

The cooling from recombination, (see e.g. Schwarz 
et al 1972), is given by: 

T)'T 

La = 

= - a:) - (6) 

mp J 
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where the constants Bi, B 2 , C, (5, A. 2 s,is are defined below. 
The function I (see Eq. 17) is the collisional ionization (see, 
e.g., Defouw 1970). Also, the term takes into account 
the photoionization and the ionization due to collisions. 

The Compton cooling-heating is given by (see, e.g., 
Peebles 1968): 


ic = 


4:kNa^ bT^x 

rrieC 


{Tm — Tr). 


(7) 


At high redshifts the Compton cooling-heating provides 
an important mechanism of energy exchange between the 
CBR and the electrons. 

For the cooling by molecular hydrogen, , we follow 
Lepp & Shull (1983), who give the following equations: 


c) H + hu ^ H + e 

C 3 = 1.5 X 10-2j;2.4g-(8.75xl03/TC g-1 (MST) 

d) H+ + R- ^ 2H + hv 

C4 = 1.7 X 10“® X cm^s-^ (PAMS) 

e) H+ + ^ H+ + e 

C5 = 5.6 X 10-®r-° ®25 (PBCMW) 

i) H + H+ ^ Hf + hv 
C6 = 3.4 X cuPs-^ (RP) 

g) Hf + H ^H2+ H+ 

C7 = 6.4 X 10"“ (KAH) 


L 


H2 


A 

P 


(erg g ^) 


h) H+ + hu ^ H+ + H 
eg = 1.1 X s"^ (MST) 


The electronic ground state of these molecules may be 
excited to rotational and vibrational levels as a result of 
collisions with atomic hydrogen, emiting infrared photons 
and hence cooling the cloud. We use the equations given 
by Lepp & Shull (1983) where they represent a good fit for 
WOK T™ s; 10®Tf. 

The equation to determine A is given by: 


i) H+ +e^ 2H+ 

C9 = 1.35 X 10“'^r"^^^ (GBD) 

i) H2+H^3H 

CIO = + (LS) 


A = n 


H2 


Lrh 


1 + 


Lrh 

Lrl 


+ ■ 


Lvf 


1 + 


Lvh 
Lvl 0 


(9) 


with: 

(i) L^„ = l.lxl0-“e-(®"""7^-); 


(ii) = 8.18 X g-(iooo/T„,) 

(for Tm > 1635A:) 

or 

= 11.45 X 
(for Tm < 1635A'); 

(iii) = 3.9 X lQ-i9^-(611S/Tm) 

(for Tm > 1087A:) 

or 

^ 2q(-19.24+0.47j/-1.247s/2) 

(for Tm < 1087A'); 

(iv) = {n°/J + 1.2n°/'^) x 1.38 x io-22g-(9234/T„) 
(for Tm ^ 4031K) 

or 

Lrl = + 1 . 2 n°„"") X 10 (- 22 . 9 -O. 553 ,- 1 . 148 j, 2 ) 

(for Tm < 4031A'). 


In the above equations is the numerical density of 
H 2 , Ugf is the numerical density of H and y — log(rr/10'*). 

The reactions which produce H 2 , used in our calcula¬ 
tions, are shown below: 


where a are the rates of production (or destruction) of the 
ions or molecules, 

= 10(4-OA16s-0.327s^) ^^3^ ^ ^ log(r^/10^), 

= 3.52 X cm^s-\ 

= 6.11 X i0~i4e-(29300/T„) g^3g-i jf 7390 a: 

or 

fcg = 2.67xl0"^®e"(®’’®°/^'")'' if < 7390^. 

The references used for the above rates are: (H) - Hi- 
rasawa 1969, (BD) - Browne & Dalgarno 1969, (MST) - 
Matsuda, Sato & Takeda 1969, (PAMS) - Peterson, Abertz, 
Moseley & Sheridan 1971, (PBCMW) - Poulaert, Brouil- 
lard, Claeys, McGowant & Wassenhove 1978, (RP) - Ra- 
maker & Peek 1976, (KAH) - Karpas, Anicich & Huntress 
1979, (GBD) - Giusti-Suzor, Bardsley & Derkits 1983, (LS) 
- Lepp & Shull 1983. 

Due to the fact that the clouds here studied are made of 
pure hydrogen, the above reactions are the most relevant to 
the formation and destruction of molecular hydrogen. Even 
if we included He, or even Li, in our studies, they would not 
contribute significantly to the cooling of the cloud, due to 
the temperatures and densities involved in our calculations. 

The rates and production of ions and molecules are cal¬ 
culated by: 

a) Molecular formation rate via H~ 


a) H + e ^ H + hv 
Cl = 1.1 X W-^^Tm (H) 


x{l — x)nci 

n[(l - x)c 2 + x(c4 + cs)] -I- C3 


( 10 ) 


b) R + R- ^R2 + e 
C2 — 2 X 10“® cm^s~^ (BD) 


dx. 


dt 


(c 2 X„_ - ciox„Jn(l - x) 


( 11 ) 
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b) Molecular formation rate via H 2 

_ C5 - (1 + x)Cfi\ 

^2 n[(l — a;)c7 + a:c9] + cs 


For the degree of ionization of the Universe, we use di¬ 
rectly the results of the study done by Peebles (1968), there¬ 
fore without taking into account the collisional ionization, 
since it is not important for the matter evolution of the Uni¬ 
verse. Thus, the degree of ionization is given by: 


= n(l - a;)[c 7 a;^+ - dox^J (13) 

The Lyman-a cooling is (e.g., Calberg 1981) given by: 

Lc = 1.25 X (erg/s) (14) 

/I 27 + O 2 I 

where 4.27 = 8.272 s“^ is the two photon emission rate, 
C 21 = 1.2 X s“^ the collisional de-exitation 

rate and C 12 — 2(7216“^^'^®^^°^ 

For the ionization degree we follow Peebles (1968) mod¬ 
ified to take into account the collisional ionization. The ion¬ 
ization rate is then written as follows: 


^ = cl-x)- ^]+I (15) 

Dt I nip j 

where Bi, is the bound energy of the ground state and B 2 
is the bound energy of the hrst excited state and: 

A2s,1s + P h'^ 

with A 2 s,is = 8.227 s“^ and a — 2.84 x cmPs~^, 

the recombination rate. 

The collisional ionization rate / (see e.g., Defouw 1970) 
is given by: 


I = 1.23 X 10“®a:fV(l -b 

Bi 


(17) 


3 BOUNDARY CONDITIONS 

The clouds evolve within a medium (the Universe) of density 
pu and thus, we take into account the influence of the ex¬ 
ternal medium on the evolution of the density perturbation. 
The expansion of the Universe is also considered because 
the physical parameters of the primordial clouds such as the 
density, the temperature and the degree of ionization de¬ 
pends on the Universe expansion rate. It is worth stressing 
that in many studies of structure formation the influence of 
the external medium on the structure under study was not 
taken into consideration. 

In order to calculate the temperature of the matter of 
the Universe we consider, in the cooling function, only the 
contribution due to the Compton cooling-heating and the 
cooling due to the recombination processes. The other phys¬ 
ical processes considered for the cloud are not relevant for 
the evolution of the matter temperature of the Universe. 

The energy equation is then written as 


dEu B dpu 

dt Pu dt 


(18) 


^ - Xu) - (19) 

at I rrip J 

^ ^ A2s,i. ^ (27rmefcr.)®/® ^-(B,/fcT.)^ (20) 

A2s,1^+P /l3 

with A 2 s,is = 8.227 and a = 2.84 x cm^s~^. 

4 CALCULATIONS AND DISCUSSIONS 

The formation of objects in the Universe may be produced 
directly either from primordial density perturbations or 
through fragmentation or still agregation of objects that are 
also produced from primordial density perturbations. 

We study here, in particular, the formation of the Pop¬ 
ulation III objects ranging from WMq to 10® M©. In the 
present work we study the formation of these objects consid¬ 
ering that they could be produced directly from a primordial 
density perturbation spectrum. 

In many studies dealing with structure formation, it is 
considered that the proto-object is already detached from 
the rest of the Universe and as a result does not undergo 
the influence of the expansion of the Universe. These calcu¬ 
lations do not take into account how a density perturbation 
becomes a coupled object, and the studies begin at a redshift 
where a couple of the physical processes are not important 
any more. 

The structures that we observe today certainly depend 
on the kind of the density perturbations produced primor- 
dially, and also on the physical processes present during and 
after the recombination era that are very important in de¬ 
velopment of these density perturbations. 

What we do here is to study one step before the kind of 
study considered elsewhere, we mean, to see how a density 
perturbation could become a proto-object. In so doing, we 
start our calculations from the recombination era, and take 
into account a series of physical processes present during and 
after the recombination era. The evolution of a given density 
perturbation depends also on, the kind of density perturba¬ 
tion spectrum and on the kind of cosmology considered (i.e., 
values of the Hubble constant and density parameters for 
baryonic and non-baryonic dark matter). 

To see how the formation of proto-objects depend on 
the particular spectrum used we consider different kinds of 
spectra. Also, we consider different cosmologies, restricting, 
however, to those cosmologies without non-baryonic dark 
matter. In another article (Oliveira et al 1998, hereafter pa¬ 
per II) we study structure formation taking into account the 
presence of non-baryonic dark matter. 

Still, to see how the physical processes contribute to 
the evolution of a given density perturbation, particularly 
those effects in general not considered by other authors, we 
performed some calculations with and without such effects 
(below we see in details which physical processes we are 
referring to). 
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Concerning the presentation of our results, we show, in 
particular, three points in the cloud: the centre (that is the 
first shell), the middle shell, and the last shell (that is the 
edge of the cloud) as a function of time. In order to analyze 
the radial behavior, we show five different times: ti - time 
at the recombination era, t 2 = (ts — ti)/2, ts = (tf — ti)/2, 
t 4 = {tf — t 3.)/2 and tf - the time at the collapse of the core. 

Also, as mentioned before, all the calculations were 
started at the beginning of the recombination era {zrec ~ 
1500). Due to the expansion of the Universe, even the clouds 
that are already Jeans unstable and eventually collapse, un¬ 
dergo, initially, an expansion phase before collapsing. It can 
be shown (see, e.g., Peebles 1993 and Coles & Lucchin 1995) 
that for a cloud which stops expanding, its mean density rel¬ 
ative to the background is given by: 



The above figure corresponds to: 


^ j ~ 4.6 (22) 

It is worth stressing that the above number is a lower 
limit since it depends only on the expansion rate of the Uni¬ 
verse (in this particular case an Einstein - de Sitter Uni¬ 
verse) and on the gravity of the cloud. For the cases in which 
other physical processes are important the clouds stop ex¬ 
panding for 

> 4.6, (23) 

far greater than the linear regime. In this way there could 
exist clouds, in the nonlinear regime, that neither collapse 
nor stop expanding (see also paper II). 

In the work by de Araujo & Opher (1988,1989) it is de¬ 
fined the collapse redshift when the cloud radius has 1/100 of 
the turn around radius (when the velocity u = 0, i.e., when 
the cloud stop expanding). In our calculation we consider 
that the cloud collapses when Sp/p ^ 10^. It is worth stress¬ 
ing that what we call a collapse it is, in fact, a condition for 
which the collapse of a proto-object become irreversible. 

Let us have a look at the models for a Universe with 
fl = 0.1 and h = 1.0 (the Hubble constant in unit of 100 
km s“^Mpc“^), and let us consider a density perturbation 
given by Eq.(l) with Mo = 10^®M©. 

Let us first have a look at the evolution of the density 
contrast (Fig.l) that gives us information about the collapse 
of the clouds. Generally speaking, we see that there is no 
collapse for M < lO'^M©. 

As expected, for lOM© the density contrast (as seen 
in Fig.l) falls rapidly in all shells and almost completely 
disappears when z ~ 1250. For M = lO^M© the same hap¬ 
pens when 2 ~ 1300. For lower redshifts the density con¬ 
trast oscillates for both masses. These results are expected, 
because at the time of the recombination era, these clouds 
were not Jeans unstable. Still, they tend to oscillate, but due 
to photon-drag the perturbations are strongly dissipated. 

For the masses 10®M© and lO^M© the fate is a little 
bit different. In these cases the perturbations oscillate for 


different shells, as can be seen, for instance, in the inter¬ 
mediate and internal shells of M = 10® M©. For the exter¬ 
nal shells the perturbation decreases almost without oscilla¬ 
tion. For the mass 10^ M© the behavior is similar to that of 
M = 10®M©. Again, due to the fact that these clouds are 
not Jeans unstable at the time of the recombination era, the 
perturbations are strongly dissipated by the photon-drag. It 
is also worth noting that the cooling mechanisms are not ef¬ 
ficient enough to cool down the clouds and turn them Jeans 
unstable. 

The masses 10® M© and 10® M© suffer partial collapse, 
i.e., the core collapses while the external shells oscillate. The 
internal shells of, for example, M = 1O®M0, initially oscil¬ 
late but, due to the cooling mechanisms, this part of the 
cloud becomes Jeans unstable and then collapses. 

It is worth stressing that we start all the calculations 
with a “top hat” density contrast profile. For the masses 
lOM© and 10^ M©, the density contrast remains almost uni¬ 
form as a function of the radial coordinate. For the masses 
10®M© and lO^M© small oscillations occur along the radial 
coordinate. For the masses 10® M© and 10® M©, the den¬ 
sity contrast profile begins approximately uniform, but it 
does not remain so as time goes on. This behavior shows 
that the “top hat” profile is not maintained throughout the 
cloud evolution. A “top hat” profile is maintained only when 
forces proportional to the spherical coordinate appears. As 
soon as the pressure gradients begin to be important in the 
cloud evolution the “top hat” density profile is destroyed. 

Alternatively in Fig.2 we illustrate the time evolution of 
the density of the cloud. We see, clearly the decrease in the 
density as time goes on for clouds with masses smaller than 
lO'^M© and the partial collapse underwent by the clouds 
with 10® M© and 10® M©. We also show the morphological 
behavior of the density profile (see Fig.3) in order to see how 
it is modified as a function of time. 

Concerning the velocity profile, it begins linear but it 
changes with time. It is worth stressing that a linear ve¬ 
locity profile is consistent with a “top hat” density profile 
and vice-versa. The masses lOM© to 10^M©, for example, 
oscillate accompanying the similar behavior of the density 
contrast. For the masses 10®M© and 10®M© the velocity 
decreases when the density contrast increases, as expected. 
We conclude that the linear profile of the velocity, that, e.g., 
de Araujo & Opher (1989) assumed, appears to be a good 
approximation for the masses lOM© to 10“^M©. But for the 
masses 10® M© and 10® M© the approximation is not valid. 
The internal pressure and the physical processes work to de¬ 
stroy the velocity linear profile, as a result the density profile 
is not a “top hat” like any more. 

In the study of the collapse, the thermal and chemi¬ 
cal history of the clouds are worth studying due to the fact 
that the evolution of the clouds in what concern, for exam¬ 
ple, their fragmentation depends on how the temperature 
of the clouds evolve which, by other hand, depend on their 
chemical ingredients. For primordial clouds, in particular, 
the thermal history will also depend on the Compton cool¬ 
ing (heating) processes, that maintains the temperature of 
the cloud close to the temperature the CBR. This cooling 
processes is particularly relevant for high values of z (i.e., 
2 > 300). 

As we explained above the heat conduction and the 
Bremsstrahlung is not taken into account. It is worth stress- 
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ing that in our study we do not take, as initial condition, a 
virialized object as considered by some authors. Our calcu¬ 
lations start at the recombination era when the temperature 
of the matter inside the perturbation has the same value of 
that of the radiation. For clouds that collapse at high red- 
shifts, in general, the principal mechanism that acts in the 
evolution and collapse of the perturbations is the Compton 
heating-cooling (see results presented in Tables 5-8). Thus, 
during the expansion phase of the perturbations, when the 
Compton heating-cooling is efficient, the temperature of the 
matter inside the perturbation is almost the same tempera¬ 
ture of the radiation. 

At the redshift collapse of the clouds, the temperature 
of the matter increases by a factor 4 or 5 in relation to 
its correspondent value at the turn-around. Thus the tem¬ 
perature of the clouds that collapse, due to shocks and the 
collapse itself, can rise to ~ \(fK depending on the power 
spectrum and on the mass of the cloud (see, e.g. Table 2). 
Also, depending on the power spectrum and on the mass of 
the cloud, the shocks and the collapse itself are not able to 
rise the temperature of the cloud to ^ lO'^A (see, e.g.. Table 
3 and paper II). 

Certainly, these results are dependent on the particular 
way in which we defined the collapse of the clouds. That 
is, we consider that an object collapsed when the density 
contrast is greater than 10^. The results also depend on the 
thermal history of the clouds (that depends on the physical 
processes here considered). As a result no strong shocks oc¬ 
cur that could eventually heat the clouds to temperatures 
> 10'*K. If we follow the calculations to higher values of 
density than the established by our collapse criterion, the 
temperature will rise above the values here obtained. It is 
possible that to density contrasts higher than 10® strong 
shocks occur and they would probably heat the gas to val¬ 
ues of the order of the virial temperature, as a result other 
physical processes, as e.g. Bremsstrahlung, would be impor¬ 
tant to the evolution of the clouds. 

Let us examine the temperature, pressure and molecular 
density formation of the clouds studied. The temperature, 
pressure and molecular density fall with time for the masses 
ranging from lOM© to IO^Mq, as seen in Fig.4. These are, 
in fact, expected results since the density perturbations are 
disappearing as time goes on. As a result, the cloud temper¬ 
ature, for example, tends to be equal to the temperature of 
the matter of the Universe. 

For clouds with masses ranging from 10 ®Mq to 10®M©, 
the behavior is very different. The pressure, for example, 
grows, since the cloud begins to collapse. The temperature 
initially decreases, but as time goes on it increases, but not 
significantly, which is so due to the fact that the cooling 
mechanisms become important and cool the cloud. During 
the collapse process there is an efficient formation of molec¬ 
ular hydrogen. The molecular hydrogen, and also the Comp¬ 
ton cooling process are the most important cooling mecha¬ 
nisms for primordial clouds. It is worth mentioning that the 
Compton cooling (heating) process is important for high red- 
shift (i.e. a > 300) as above mentioned, while the molecular 
hydrogen cooling is important for 2 < 300, which is so due 
to the fact that the CBR photoionize the H~ and H 2 ions 
that intermediate the formation of molecular hydrogen. 

An analysis of the molecular density versus tempera¬ 
ture of the cloud is shown in Fig.5. For the cloud of 10®M©, 


for example, the molecular density increases when the tem¬ 
perature is in the interval lOOA <T < SOOA, which occurs 
during the cloud collapse. The molecular hydrogen is very 
important in the cooling processes of the cloud, as already 
mentioned. 

For the masses lOM© to 10"^M© the pressure proHle 
remains constant with time along the radial coordinate; 
only tiny variations are seen for the mass lO^M©. For other 
masses, the profile changes rapidly with time along the ra¬ 
dial coordinate. The production of the molecules, and the 
consequent cooling of the cloud, is always more important 
in the central part of the cloud. 

In Fig.6 we show how the radii of the shells of the clouds 
change with time. For the masses 1O®M0 and 10®M©, for 
example, we clearly see how the different shells of the clouds 
evolve with time. It is also possible to see that these clouds 
undergo a partial collapse, the outer most shells do not stop 
expanding. 

In Table 1 we present, in particular, in its second and 
third column, a comparison between our results and the re¬ 
sults of de Araujo (1990) for Mo = 10'^®M©. There we ob¬ 
serve, e.g., that the mass 10®M© has collapse redshift simi¬ 
lar to that obtained by de Araujo (1990). The mass 10"^M©, 
however, collapses at 2c = 5.4 in the models of de Araujo 
(1990) but here there is no collapse for such a mass. The 
mass 10®M© collapses at 2c = 260 according to de Araujo 
(1990), but we obtained a collapse at Zc = 22. For clouds 
with masses larger than 10® M© the collapse redshift is al¬ 
most the same, for these cases the pressure is not very im¬ 
portant in the evolution of the cloud, this explain why the 
results are the same (see also paper II). For those cases in 
which the internal pressure is important on the evolution of 
the cloud the results are too different. In conclusion, it is not 
a good assumption to maintain a “top hat” density profile 
for clouds with masses less than 10®M©, as de Araujo (1990) 
considered throughout their calculations. 

As mentioned earlier the clouds stop expanding satisfy¬ 
ing Eq.(23), therefore well above the linear regime. Again we 
see in Table 1 in particular for a cloud of 10® M© that the 
present results agree with those by de Araujo (1990). For 
clouds in which the pressure is important the turn around 
occurs for values above the value given by Eq.(23). In the 
paper II we discuss this issue for clouds with masses larger 
than 10® M©, where we also include the presence of non- 
baryonic dark matter. 

Another interesting comparison has to do with the min¬ 
imum mass that can collapse in our model. We have found 
that the minimum mass is ~ 8.5 x 10^M© (the collapse red¬ 
shift is 12), that is almost one order of magnitude larger 
than that found by de Araujo (1990). This shows again that 
to fix the density profile is not a good assumption. 

We obtained also a mass limit that does not collapse 
but remain bounded, namely, 3.8 x 10^M©. This bounded 
cloud reaches a density contrast of 10^, but the collapse does 
not evolve to a denser stage, it oscillates until the present 
time. In Fig.7 we see the behavior of the radii of the shells 
as a function of time for this mass. 

It is worth stressing that these results just discussed 
are for a Universe model with h = 1 and — 0.1, with a 
reference mass Mo = 10^®M©. 

Obviously the results presented here depend on the par¬ 
ticular spectrum of perturbations used. The minimum mass. 
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Table 1. The redshifts of collapse. 


M/M© 

2c 

2^ 

^c 

z* 

°ta 

Sta 

10 ^ 

NC 

5.4 

NC 

8.0 

— 

4.4 x 10^ 

- 

- 

40 

- 

- 

8.5 x lO"' 

12 

- 

498 

- 

6.3 

10 ® 

22 

260 

477 

6.7 

5.7 

10 ® 

197 

180 

643 

5.1 

5.0 


The redshifts of collapse of the present paper, z^, and 
de Araujo’s studies (de Araujo 1990), zf, for Mo = 

10 ^® Mq; and redshifts of collapse of the present paper 
for Mo = 4 X 10^'^Mq, z*. The density contrast at the 
turn around for de Araujo’s studies, and for the 
present paper, &ta, for Mo = NC stands for 

“no collapse”. 

for example, that collapses would be different if we used a 
different spectrum of density perturbations. In Table 1 we 
also include the collapse redshift for the spectrum normal¬ 
ized with a reference mass Mo —Ax (instead of 

The density perturbations with Mo = 4 x 
are seven times greater than those with Mo — IC^Mq. 
We may see that the collapse redshifts are altered signifi¬ 
cantly, and also the minimum mass that collapses is now 
4.4 X ICMq with a collapse redshift of ~ 40 (see Table 2). 

Due to the fact that the collapses occur at a higher 
redshift for a large Mo, there is also a modification in the 
importance of the physical processes in the evolution of the 
clouds. The molecular hydrogen cooling, for example, is less 
significant for collapses which occur at high redshift. 

When we change the content of baryonic matter from 
Q, — 0.1 to = 0.2 and maintain h — 1 and Mo = IO^^Mq 
the minimum mass becomes 5 x lO^M©; the collapse occur- 
ing at z=9 (see Table 3). In Table 3 we show the collapse 
redshifts for clouds of different masses. 

For h = 0.5, 14 = 0.1 and Mo = lO^^M© we obtain that 
2 X 10® M© is the minimum mass that collapses. The collapse 
occurs at 2=0.78 (see Table 4). 

The above results show that the minimum mass does 
not depend strongly neither on the kind of spectra used 
nor on the values of the Hubble constant and the density 
parameter assumed. By other hand, the collapse redshifts 
depend strongly on the Universe model and on the spectra 
adopted. 

As we have already mentioned, the present study in¬ 
cludes a series of physical processes not in general considered 
in the literature, namely, the Compton cooling (heating), the 
photon-drag due to the cosmic background radiation are ex¬ 
amples. 

We thus studied the influence of the cooling-heating 
processes and the photon-drag on the evolution of the clouds 
here considered. We chose, in particular, for this study, the 
clouds of 3.8 X 10^, 8.5 x 10^, 10® and 10®M©, the results are 
presented in Tables 5-8 respectively. To do this we studied 
individually the effects of all the four cooling-heating pro¬ 
cesses and the photon drag. The values in these tables have 
been taken for z^, that is, the redshift when the collapse 
occurs. 

First let us analyze the cooling from recombination. For 
the mass 3.8 x lO^M© we again do not have collapse when 
we disregard this cooling mechanism; for 8.5 x lO^M© the 
collapse is delayed. The collapse redshift changes from 2 ~ 


Table 2. Values of the variables 
£7 = 0.1 and Mo = 4 x 1O^^M0. 

when we 

used h = 1 . 0 . 

M/M© 


^H2 Cm ®) 

T{K) 

p {g cm~^) 

lO"' 

NC 

_ 

— 

— 

4.4 X 10"^ 

40 

3 X 10-2 

7 X 10® 

8 X 10“^® 

10 ® 

48 

3 X 10-2 

7 X 10® 

8 X 10“^® 

10 ® 

624 

2 X 10-2 

7 X 10® 

5 X lO-^"^ 


Table 3. Values of the variables 
n = 0.2 and Mo = 1 O®®M 0 . 

when we 

used h = 1 . 0 . 

M/M© 

2^0 

Cm ®) 

TiK) 

p {g cm-^) 

5 X 10^ 

9 

300 

100 

3 X 10-22 

8.5 X 10^ 

90 

100 

500 

3 X 10-1® 

10® 

123 

320 

500 

1 X 10-1® 

10® 

207 

500 

10® 

3 X 10-1® 


Table 4. Values of the variables when we used h = 0.5, 
n = 0.1 and Mo = 10 ^^Mq. 


M/M© 


{cm ®) 

T{K) 

P 

(9 

cm 

-3 

8.5 X 10^ 

NC 

2 X 10-12 

10 

3 

X 

10 “ 

■31 

10 ® 

NC 

3 X 10-1® 

10 

3 

X 

10 - 

■31 

2 X 10® 

0.8 

3 X 10-® 

50 

3 

X 

10 “ 

■25 

10 ® 

158 

63 

320 

2 

X 

10 - 

■19 


12 with all processes to 2 ~ 6.8 without the recombination 
cooling. For the masses 10® M© and 10® M© this cooling does 
not affect the collapse redshift. 

The Compton cooling-heating works against the col¬ 
lapse for the masses from 3.8 x 10^M© to 10®M©. When 
a cloud (density perturbation) begins to evolve, at the be¬ 
ginning of the recombination era, it first expands and then 
stop expanding if it becomes Jeans unstable, as already men¬ 
tioned. During the expansion phase it would cool down adi- 
abatically if we did not have any cooling-heating processes. 
In this way, the cloud could stop expanding and collapses 
earlier. Due to the Compton heating the cloud is heated to 
temperatures near to that of the cosmic background radi¬ 
ation, as a result the cloud stop expanding and collapses 
later. 

Also, disregarding the Compton heating a cloud of 3.8 x 
lO^M© collapses. This means that the minimum mass de¬ 
pends on this very processes. For masses with M < 10®M© 
the collapse occurs earlier as we may see in Tables 5-8. For 
clouds with M > 10®M©, however, the collapse redshift is 
not affected, again due to the fact that the pressure for this 
cloud does not significantly affect the collapse, because these 
clouds become Jeans unstable just after the recombination 
era. 

In Fig.8 we present the time evolution of the tem¬ 
perature, pressure and molecular hydrogen density, with 
all physical processes included and without the Compton 
cooling-heating processes. It is seen that the evolution of 
the above physical quantities are very different, even for the 
case in which the redshift collapse is not significantly al¬ 
tered, namely, for clouds of 10®M©. In conclusion, the ther¬ 
mal history of the clouds studied is strongly dependent on 
the Compton heating process. 
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Table 5. Values of the 

variables for 

M = 3.8 X 

10^ Mq with 

h = 1.0, 

= 0.1 and 

Mo = 10*® M©. 








with all 

without 

without 

without 

without 

without 

Variable 

processes 

Ln 

Lc 


Let 

photon drag 


NC 

NC 

156 

NC 

NC 

NC 


10“® 

10“® 

102 

- 

10-® 

10-10 

T (K) 

50 

50 

102 

50 

50 

50 

P {g crn ®) 

3 X 10-27 

3 X 10-27 

3 X 10-1® 

2 X 10-27 

3 X 10-27 

1 X 10-®° 


Table 6 Values of the 

variables for 

M = 8.5 X 

IO^Mq with 

h = 1.0, 

= 0.1 and 

Mo = 10*® AJ©. 








with all 

without 

without 

without 

without 

without 

Variable 

processes 

La 

Lc 


Let 

photon drag 


12 

6.8 

252 

2.7 

6.8 

2 


10-2 

10-2 

3 X 10-2 

- 

10-2 

3 X 10-* 

T (K) 

102 

102 

7 X 10® 

4 X 102 

102 

102 

p {g cm, ®) 

1 X 10-22 

1 X 10-22 

6 X 10-2° 

1 X 10-2® 

1 X 10-22 

3 X 10-24 


Table 7. Values of the variables for M = 10® Mq with h = 1.0, Q = 0.1 and Mo = 10^® Mq. 


Variable 

with all 

processes 

without 

La 

without 

Lc 

without 

without 

Let 

without 
photon drag 


24 

24 

220 

5 

24 

13 


0.3 

0.5 

102 

- 

0.3 

6 X 10-2 

T (K) 

200 

200 

2 X 10® 

200 

200 

102 

P [g cm ®) 

3 X 10-21 

3 X 10-21 

2 X 10-1® 

3 X 10-21 

3 X 10-21 

6 X 10-22 


Table 8 Values of the variables for M 

= with h = 1.0, £7 

= 0.1 and Mo 

= 1O1®M0. 


with all 

without 

without 

without 

without 

without 

Variable 

processes 

La 

La 


Let 

photon drag 

2^ 

194 

194 

196 

157 

194 

246 

Tiffs (cm"®) 

800 

800 

800 

- 

102 

102 

T {K) 

800 

800 

10® 

8 X 10® 

10® 

1300 

p {g cm 

2 X 10-1® 

2 X 10-1® 

2 X 10-1® 

1 X 10-1® 

2 X 10-1® 

3 X 10-1® 


If we disregard the cooling by molecular hydrogen, the 
mass 3.8 x IO^Mq is not signihcantly affected, there is no 
collapse in this case, as already seen there is no collapse 
even considering the molecular hydrogen cooling. The col¬ 
lapse for the other masses occur for lower values of redshift. 
Without such an important cooling process the pressure is 
greater and as a result there is a delay in the collapse of the 
cloud. In Fig.9 we can see the different behavior of the tem¬ 
perature, pressure and molecular density when this cooling 
is disregarded. 

For the Lyman-a cooling the behavior of the collapse is 
similar to the cooling due to the recombination. 

If we disregard the photon-drag process we observed 
that the collapse is delayed for masses 8.5 x 10'*M© and 
10®M© and the collapse occurs earlier for the mass 10®M©. 

The photon drag works like a brake for the expansion 
and for the collapse. If we have masses around the Jeans 
mass the behavior of these masses near the epoch of the 
collapse is strongly influenced by the photon drag. 

When we disregard the photon drag for masses smaller 


than 10® M© the collapse is delayed, due to the fact that the 
cloud stop expanding latter, beginning as a result to collapse 
latter. This occurs, also because these clouds have M ^ Mj 
(the Jeans mass). 

If we have masses larger than 10®M© (Jeans unstable 
masses) when we disconsider the photon drag, the cloud 
stop expanding earlier, and due to the fact that the pressure 
is not so important for these cases, the cloud collapses in 
a. t ^ tf[ (free-fall time). As a result the collapses occur 
earlier, without such an effect. 

As a general conclusion we see that the cooling-heating 
processes are less effective for the mass 10® M©, in partic¬ 
ular, concerning the epoch of the collapse, due to the fact 
that the pressure of the cloud in this case does not affect sig¬ 
nificantly its evolution. By other hand, the thermal history 
of the cloud, as well as, the amount of molecular hydrogen 
produced depend signihcantly on the physical processes here 
considered. 

An important point concerning the present study is that 
we take into account the appropriate physics and begin the 
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calculations in the appropriate era, i.e., at the recombination 
era. These facts strongly influence the evolution of the clouds 
as we have shown above. 

5 CONCLUSIONS 

In the present work we studied the evolution of density per¬ 
turbations with masses within the range 10 — 10® Mq taking 
into account a series of physical processes which are present 
during and after the recombination era. In order to perform 
such a study a spherical Lagrangian hydrodynamical code 
was written, with which it was possible to follow the spatial 
and time evolution of the density perturbations. 

Our main conclusions are: 

a) For clouds with M < IC'Mq, no collapse occurs. In 
this case, the density perturbations are strongly dissipated 
by the photon-drag. Only a residual density perturbation 
remains; this conclusion holds for the two density perturba¬ 
tions spectra and for the different combinations of h and 
adopted; 

b) Clouds of 10® and 10 ®Mq present a partial collapse; 

c) The minimal mass which collapses does not change sig¬ 
nificantly when we alter h, Q and Mq, namely, Mmin ~ 
10^ - 10® M©; 

d) When we disregard, for example, the Compton cooling¬ 
heating, the collapse for masses > 8.5 x ICMq occurs ear¬ 
lier; 

e) The photon-drag is a very important process, the col¬ 
lapse redshift, for example, is changed significantly when we 
disregard it, but it works differently for higher or smaller 
masses. For the smaller masses the photon-drag dissipates 
their density perturbations and for the higher masses it 
causes a delay in the collapse of the clouds. 

ACKNOWLEDGEMENTS 

The calculations were performed on a HP-Apollo 9000 (pur¬ 
chased by the Brazilian agency FAPESP) and on the CRAY 
computers EL98 and J90 (CCE/USP). 

We would like to thank Dr David Weinberg, the referee, 
for his helpful suggestions and comments that greatly im¬ 
proved the present version of our paper. We would also like 
to thank the Brazilian agencies CAPES (SRO) and CNPq 
(ODM, JCNA and RO) for financial support. 

REFERENCES 

Browne J.C., Dalgarno A., 1969, J. Phys, 2, 885 
Calberg R.G., 1981, MNRAS, 197, 1021 
Cen R., 1992, ApJ, 78, 341 

Coles P., Lucchin F., 1995, Cosmology. The Origin and Evolution 
of Cosmic Structure, John Wiley & Sons, New York 
de Araujo J.C.N., 1990, PhD thesis, ACA-120, lAC-USP 
de Araujo J.C.N., Opher R., 1988, MNRAS, 231, 923 
de Araujo J.C.N., Opher R., 1989, MNRAS, 239, 371 
de Araujo J.C.N., Opher R., 1994, ApJ, 437, 556 
Defouw R.J., 1970, ApJ, 161, 55 

Giusti-Suzor A., Bardsley J.N., Derkits C., 1983, Phys. Rev. A, 
28, 682 

Gott III R., Rees M.J., 1975, A&A, 45, 365 
Haiman Z., Thoul A.A., Loeb A., 1996, ApJ, 464, 523 
Hirasawa T., 1969, Prog. Theor. Phys., 42, 523 


Hutchins J.B., 1976, ApJ, 205, 103 

Karpas Z., Anicich V., Huntress W.T., 1979, J. Chem. Phys., 70, 
2877 

Lepp S., Shull J.M., 1983, ApJ, 270, 578 

Matsuda T., Sato H., Takeda H., 1969, Prog. Theor. Phys., 42, 
219 

Miranda O.D., 1998, in preparation 

Oliveira S.R., Miranda O.D., de Araujo J.C.N., Opher R., 1998, 
MNRAS (submitted) - Paper H 

Palla F., Salpeter E., Stabler S., 1983, ApJ, 271, 632 

Peebles P.J.E., 1968, ApJ, 153, 1 

Peebles P.J.E., 1993, Principles of Physical Cosmology, Princeton 
UP, Princeton 

Peterson J.R., Abertz W.H., Moseley J.T., Sheridan J.R. 1971, 
Phys. Rev. A, 3, 1651 

Poulaert G., Brouillard F., Claeys W., McGowant J. W., Wassen- 
hove G. van, 1978, J. Phys. B, 11, 21 

Ramaker D., Peek J., 1976, Phys. Rev. A, 13, 58 

Richtmeyer R.D., Morton K.W., 1967, Difference Methods for 
Initial Value Problems, Interscience Pub., New York 

Schwartz J., MacCray R., Stein R.F., 1972, ApJ, 175, 673 

Thoul A.A., Weinberg D., 1995, ApJ, 442, 480 


FIGURE CAPTIONS 


Fig.l — The evolution of the density contrast (5) with 
time. The solid line is the internal shell, the short dashed 
line is the middle shell and the long dashed line is the 
external shell. Time is in years (the input parameters are 
Mo = 10^®M©, h = 1.0, O = 0.1). 

Fig.2 — The evolution of the mass density p {g cm ®) 
with time. The captions and input parameters for the curves 
are the same as in Fig.l. 

Fig.3 — The mass density versus the radial coordi¬ 
nate (in pc) of the cloud. The solid line is the initial time 
(ti), the short dashed line, long dashed line and very long 
dashed line are intermediate times (respectively t 2 , ts and 
t 4 ). The dashed-pointed line is the final time (t/). The input 
parameters are the same as in Fig.l. 

Fig.4 — The evolution of the temperature (T in K), 
pressure (P in dyn cm~^) and density of molecules {nH 2 in 
molecules cm~^) with time. The captions and input param¬ 
eters for the curves are the same as in Fig.l. 

Fig.5 — The evolution of the density of molecules (riH^) 
versus temperature. The captions and input parameters for 
the curves are the same as in Fig.l. 

Fig.6 — The evolution of the radial coordinate of the 
clouds with time. The captions and input parameters for the 
curves are the same as in Fig.l. 

Fig.7 — The evolution of the radial coordinate of the 
mass M = 3.8 x 10^M© with time. The captions and input 
parameters for the curves are the same as in Fig.l. 

Fig.8 — The evolution of the temperature (T), pressure 
(P) and density of molecules (riH^) with time for masses 
M = 3.8xlO'‘M©, M = 8.5xlO‘‘M©, M = 10®M© and M = 
10®M© with all physical processes included and when we 
disregard the Compton cooling-heating. The captions and 
input parameters for the curves are the same as in Fig.l. 

Fig.9 - The evolution of the temperature (T), pressure 
(P) and density of molecules {nH 2 ) for masses M = 3.8 x 
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10*Mq, M = 8.5 X W^Mq, M = 10®Mq and M = 10®M© 
with all physical processes included and when we disregard 
the H 2 cooling. The captions and input parameters for the 
curves are the same as in Fig.l. 


This paper has been produced using the Royal Astronomical 
Society/Blackwell Science TliiX macros. 
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